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CLAIMS 

1 . A method for the accurate determination of van der Waals parameters for high-precision 
determination of crystal structures and/or energies, comprising the steps of: 

- numerically simulating at least one crystal structure based on density functional 
theory (DFT) calculations combined with a potential energy term representing Van der 
Waals interactions; 

- providing reference data containing accurate information about said at least one 
crystal structure; 

- defining a deviation function (F) quantifying a deviation between said reference data 
and said at least one simulated crystal structure; 

- fitting at least one parameter of said van der Waals potential term in such a way as to 
minimize said deviation function (F); and 

obtaining the accurate van der Waals parameters from the best fit. 

2. A method according to claim 1, characterized in that said van der Waals potential term is 
defined as: 
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wherein fAjjOXB) is a damping function and the sum runs over all pairs of interacting atoms, 
and that said fitting step comprises fitting said damping function. 

3. A method according to claim 2, characterized in that said damping function is defined as 



and that said fitting step comprises fitting the parameter ta,b and/or the parameter n and/or the 
parameter c. 
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4. A method according to claim 2 or 3, characterized in that said fitting step furthermore 
comprises fitting said coefficient C^a,b. 
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5. A method according to any of the preceding claims, characterized in that said reference 
data are theoretical data obtained by Hartree-Fock calculations or Quantum Monte Carlo 
simulations. 

6. A method according to any of the preceding claims, characterized in that said reference 
data are experimental low-temperature crystal structure data. 

7. A method according to claim 6, characterized in that said crystal structure data are 
obtained by X-Ray or neutron scattering. 

8. A method for the accurate determination of crystal structures and/or energies, comprising 
the steps of: 

- providing a rongh estimate model of at least one crystal structure; 

- numerically simulating said at least one crystal structure based on density functional 
theory (DFT) calculations combined with a potential energy term representing Van der 
Waals interactions; and 

- obtaining said at least one crystal structure and/or its energy as a result of said 
numerical simulation, 

characterized in that said van der Waals potential term is obtained by the method according to 
any of claims 1 to 7. 

9. A method according to claim 8, characterized in that a plurality of polymorphic crystal 
structures are determined and ranked according to their respective energies. 

1 0. A method for the efficient numerical optimization of a molecular crystal structure using an 
advantageous crystal coordinate system, comprising the steps of: 

- providing a starting crystal lattice described by an initial coordinate system 
comprising lattice parameters and atomic positions in said crystal; 

- defining a so-called natural coordinate system and representing said starting crystal 
lattice in said natural coordinate system, said natural coordinate system comprising: 

• first coordinates describing symmetry-allowed lattice changes and 
defined in such a way that changes of said first coordinates do not cause 
changes of the molecular geometry or a rotation of molecules with respect to 
each other and leave fractional coordinates of molecular centres constant; 
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• second coordinates describing symmetry-allowed translations of said 
molecules in said crystal; 

• third coordinates describing symmetry-allowed rotations of said 
molecules in said crystal; 

• fourth coordinates describing symmetry-allowed changes of the 
molecular geometry; 

- transforming coordinates from said natural coordinate system to said initial 
coordinate system; 

" calculating the lattice energy and energy derivatives with respect to said initial 
coordinate system; and 

- transforming said energy derivatives from said initial coordinate system to said 
natural coordinate system, 

wherein a minimization algorithm is used for minimizing said lattice energy with respect to 
said natural coordinate system. 

11. A method for the energy ranking of polymorphic crystal structures, comprising the steps of: 

- providing rough estimate models of each of said crystal structures; 

- numerically simulating each of said crystal structures based on density functional 
theory (DFT) calculations combined with a potential energy term representing Van der 
Waals interactions. 

- obtaining accurate crystal structures and energies as a result of said numerical 
simulation; and 

- ranking said accurate crystal structures according to their respective acccurate 
energies. 

12. The method according to claim 11, characterized in that the crystals are crystals of 
pharmaceutical compounds. 

13. The method according to claim 12, characterized in that it is applied to identify the most 
stable polymorphic form of a pharmaceutical compound. 

14. The method according to any of claims 1 1 to 13, characterized in that it uses a method for 
the efficient optimization of the molecular crystal structure according to claim 10. 
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15. A computer program product comprising computer readable code for enabling a computer 
to perform a method according to any of the preceding claims when said code is executed. 
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